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^ Abstract 

^ We give new and explicitly computable examples of Gibbs-non- 

i ^ i Gibbs transitions of mean- field type, using the large deviation ap- 

proach introduced in [JJ. These examples include Brownian motion 
with small variance and related diffusion processes, such as the Ornstein- 
^r) Uhlenbeck process, as well as birth and death processes. We show for 

a large class of initial measures and diffusive dynamics both short- 
time conservation of Gibbsianness and dynamical Gibbs-non-Gibbs 

• transitions. 
CN 
O 



1 Introduction 

Starting from [3] dynamical Gibbs-non-Gibbs transitions have been consid- 
ered by several authors, see e.g. [2J, [ID], [S]- In these studies, one considers 
lattice spin systems started from a Gibbs measure \i at time zero and evolves 
it according to a Markovian dynamics (e.g. Glauber dynamics) with station- 
ary Gibbs measure v ^ \i. The question is then whether fi t , the time-evolved 
measure at time t > is a Gibbs measure. Typically this is the case for short 
times, whereas for longer times, there can be transitions from Gibbs to non- 
Gibbs (loss) and back from non-Gibbs to Gibbs (recovery). The notion of a 
"bad configuration", i.e., a point of essential discontinuity of the conditional 



1 



probabilities of the measure \xt is crucial here. Such a configuration r] spec is 
typically identified by looking at the joint distribution of the system at time 
and at time t. If conditioned on r) spec the system at time zero has a phase 
transition, then typically rj spec is a bad configuration. 

In the context of mean-field models, the authors in [9j started with an 
analysis of the most probable trajectories (in the sense of large deviations) 
of a system conditioned to arrive at time T at a given configuration. The 
setting of [S] is the Curie- Weiss model subjected to a spin-flip dynamics. A 
Gibbs-non-Gibbs transition is in this context rephrased as a phenomenon of 
"competing histories", i.e., for special terminal conditions x spec and times T 
not too small, multiple trajectories can minimize the rate function, and these 
trajectories can be selected by suitably approximating x spec . These special 
conditionings leading to multiple histories are the analogue of "bad configu- 
rations" (essential points of discontinuity of conditional probabilities of the 
measure at time t) in the Gibbs-non-Gibbs transition scenario. This "tra- 
jectory large deviation approach" has then been studied in more generality, 
including the lattice case, in [I]. 

In this paper, we apply the trajectory large deviation approach in several 
examples, both for diffusion processes and for birth and death processes. This 
leads to new and explicitly computable Gibbs-non-Gibbs transitions of mean 
field type. For processes of diffusion type, we first treat an explicit example 
for the rate function of the initial measure, and as dynamics Brownian motion 
with small variance or Ornstein-Uhlenbeck process. In all cases, we obtain the 
explicit form of the conditioned trajectories, and explicit formulas for the bad 
configuration and the time at which it becomes bad. In the case of general 
Markovian diffusion processes in a symmetric potential landscape, we show 
under reasonable conditions short-time Gibbsianness as well as appearance of 
bad configurations at large times. Next, we treat the case of continuous-time 
random walk with small increments, as arises e.g. naturally in the context of 
(properly rescaled) population dynamics. In that case, the Euler-Lagrange 
trajectories can be explicitly computed for some particular choices of the 
"birth and death" rates. Constant birth and death rates are the analogue 
of the Brownian motion case, whereas linear birth and death rates are the 
analogue of the Ornstein-Uhlenbeck process, but in that case the cost of 
optimal trajectories becomes a much more complicated expression. 

Our paper is organized as follows. In section[2]we introduce some elements 
or the Feng-Kurtz formalism, and define the notion of bad configurations in 
the present setting. In section [3] we treat diffusion processes with small 



variance, with an explicit form for the initial rate function. In section 3.3 



we treat the case of Brownian motion dynamics with different cases for the 
rate function of the initial measure. Finally, in section |5j we treat one- 
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dimensional random walks with small increments, such as rescaled birth and 
death processes. 



2 The Feng-Kurtz scheme, Euler- Lagrange tra- 
jectories, bad configurations 

We study Markov processes {X™ : < t < T} taking values in IR d , parametrized 
by a natural number n. This parameter tunes the "amount of noise" in the 
process, i.e., as n — > oo, the process becomes deterministic, and the measure 
on trajectories satisfies the large deviation principle with rate n and with a 
rate function of the form 

J{i)= [ L( ls ,%)ds (1) 
Jo 

This means more precisely that 

P ({X r t l : < t < T} « 7) « exp (-71^(7)) (2) 

to be interpreted in the usual sense of the large deviation principle with a 
suitable topology on the set of trajectories. The form ([I]) naturally follows 
from the Markov property. 

Notice that the form of the rate function does not depend on the choice 
of this topology. So one usually starts with the weakest topology, i.e., the 
product topology, and then, if possible, strengthens the topology by showing 
exponential tightness. See [1] for an illustration of this strategy in the context 
of theorems like Mogulskii's theorem. 

Since in this paper we are only interested in finding out optimal trajec- 
tories, i.e., minimizers of the rate function over a set of trajectories with 
prescribed terminal condition and open-end condition, we will not have to 
worry about the strongest topology in which the large deviation principle 
g) holds, but we are rather after (as explicit as possible) solutions of Euler- 
Lagrange problems associated to the rate function. 

In [S] a scheme is given to compute the "Lagrangian" L, see also [I] for 
an illustration of this scheme in the large-deviation view on Gibbs-non-Gibbs 
transitions. First one computes the "Hamiltonian" 

Jf?(p,x) = lim I e -»<P.*>j2> ne "<P.*> (3) 

n— >oo n 

where S£ n is the generator of the process {X™ : < t < T} (working on the re- 
variable, where p e R d is the "momentum" and where < ., . > denotes inner 
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product. . Under regularity conditions on J4?(p,x) (e.g. strict convexity), 
the associated Lagrangian is then given by the Legendre transform 



L(x,v) = sup (< v,p > —J^(x,p)) , (4) 
As an example, consider 



with generator 



then we have 



X, n = n- 1/2 B t 



S? n = —A 
2n 



v 2 



and associated Lagrangian 
which produces the rate function of the well-known Schilder's theorem 



L(x, v) 
k > j 2 



P(W : < t < T} « 7) ~ exp (~ £ ^ ds ) 

To proceed, we also want the initial point of our process to have some fluctu- 
ations. More precisely, we need for the starting point of our process an initial 
measure \i n (depending on n) on IR d , satisfying the large deviation principle 
with rate n and rate function i(x), i.e., in the sense of large deviations, we 
assume 

P (X n E A) = ii n (A) w exp(-n inf i(x)) (5) 

We call the triple ({X f n : < t < T},L,i) a stochastic system with small 
noise. 

We continue now with the definition of a bad configuration in this frame- 
work. This is motivated by the definition of a bad configuration in the context 
of mean- field models [S] , and can be viewed as the large-deviation rephrasing 
of "a phase transition at time zero conditioned on a special configuration at 
time T". 

Definition 2.1. Let ({X™ : < t < T},L,i) be a stochastic system with 
small noise. We say that a point b £M. d is bad at time T if the following two 
conditions hold. 

1. Conditional on X^, = b, X£ does not converge (as n — )■ 00 ) to a point- 
mass in distribution. 
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2. There exist two sequences b k — >• b, b k — >■ b and 5 > such that the 
variational distance between the distribution //(0,T; of Xq \X% = b k 
and the distribution fi(0,T;b k ) of Xq\X£ = b k is at least S for k large 
enough. 

The simplest example which follows also the most common scenario is 
where the distribution of Xq \X% = b converges to \{o~- a + o~ a ) and for c > b 
Xq \Xj, = c converges to S a / c -\ where a(c) — > a as c 4- b, whereas for c < b 
Xq \Xj, = c converges to 5 a /t c -j where a'(c) — > —a as c t b. This means that 
conditioned to be at time T at location b, the process has two "favourite" 
intial spots, which can be "selected" by approaching b from the right or from 
the left. 

This is the analogue of a phase transition, where the phases can be se- 
lected by appropriately approximating the bad configuration, see [3]. 

3 Diffusion processes with small variance con- 
ditioned on the future 

In this section we present examples where X™ is a diffusion process. We show 
also how from the large deviation approach we gain a new understanding of 
"short-time Gibbsianness" for a general class of drifts of the diffusion, or 
initial rate functions. 



3.1 Brownian motion 

To start with, we consider Brownian motion with small variance - starting 
from an initial distribution satisfying the large deviation principle (with rate 
n) with a non-convex rate function having two mimina at locations —a, a, 
with a > 0. More precisely, we consider the process 

*T = 4=^ (6) 

starting from an initial distribution fi n such that, informally written, 

P (X n G dx) = fi n (dx) w e- ni(x) dx (7) 

For i we make the explicit choice: 

i{x) = ( x 2 _ a 2 ) 2 (8) 

i.e., a non-convex function non-negative with zeros at —a, a and maximum 
at x = (i(x) with a = 2 is plotted in Figure [T]). 



5 



i 25 




\ 20 




\ 




\ / 10 




\ / 5 





-3 -2-1 12 



Figure 1: = (x 2 — a 2 ) 2 with a = 2 

This specific choice is for the sake of explicit analytic computability but 
many results are true for a general class of rate functions that have a similar 
graph with two zeros located at —a, a and a maximum at zero. 

More formally, we require that the sequence of initial probability measures 
{fin,n G N} satisfies the large deviation principle with rate function i given 
by ([8]). Such rate functions arise naturally in the context of mean-field models 
with continuous spins and spin-Hamiltonian depending on the magnetization. 

We are then interested in the most probable trajectory 7 with initial 
point distributed according to \x n and final point jt = 0. More precisely, by 
application of Schilder's theorem, the trajectory {X™ : < t < T} satisfies 
the LDP with rate function 

*W = lf isds + i(lo) (9) 
The optimal trajectory we are looking for is hence 

argmin{j^(7) : j(T) = 0} 

The Euler-Lagrange trajectories (extrema of the cost | J Q T 7 2 g?s corresponding 
to 0^(7)) are linear in t: 

>y t = A + Bt 
By the terminal condition B = —A/T. 
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The cost 0^(7) of this trajectory can then be rewritten as a function of 
the starting point 70 = A: 

£ ,t( A ) := J?{i) = A 4 -2a 2 A 2 W+l(-A/T) 2 T = A 4 +a(a,T)A 2 W (10) 
with 



a(a, T) = 




The behavior of this cost depends on the sign of a. If a > 0, then there is a 
unique minimum at A = 0, this case corresponds to 

If a < then there are two mimima A = A± given by 

A± = ±y/-a(a,T)/2 = ±^a 2 - (AT)- 1 (12) 

We thus conclude that, as n — > 00, the starting point is most probably 
for small T and most (and equally) probably A± for large T, which converges 
to ±a when T — > 00. Hence we have non-uniqueness of histories. 

Let us denote T, 0) the distribution of Xq conditioned on X% = 0. 
Then we have 

1. Small times, unique history. If T < T crii then 

lim /i(n, T, 0) = 5o- 

n— >oo 



2. Large times, non-unique history. If T > T crit then 

lim fi(n,T,0) = l,(S A + +5 A -). 

n— >oo z 



3. Limit of large times 

lim lim /i(n, T, 0) ->• hs a + <L a ). 

T— >cc n— >oo z 



Let us now condition on X£ = b 7^ 0. Then the most probable trajectory 
is still a straight line 7^ = A+Bt but now with terminal condition A+BT = b, 
i.e., -B = (6— A) /T. It has cost expressed in terms of the starting point 70 = A 

*w(A) = A 4 + a(a, T)A 2 + a 4 - h -A + %- (13) 
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This is the cost function <£o.t(^4) of (10) plus a linear term —hA + Min- 
imization of Sb t T{A) leads to the equation 

AA 3 + 2aA = ^ (14) 

We then have two cases: 



1. a > 0, i.e., T < T crit . Equation (14) has a unique real solution, corre 



sponding to a unique minimum Aj, of Ej,(A). This minimum converges 
to zero as b — > 0. Hence, is good for T < T crit . 



2. a < 0. Equation (14) has three real solutions. For b > we have 



one positive and two negative solutions. The positive solution denoted 
A(+, b, T) > a/— a/2 gives the minimum. The negative solutions corre- 
spond to a maximum and a local minimum. For b < the situation ex- 
actly the opposite: the unique negative solution A(—, b, T) < — a/— a/2 
correspond to the global minimum whereas the two positive solutions 
give a maximum and a local minimum. Hence is bad for all T > T crit 

In particular, for the T — > oo the positive, resp. negative minimum of the 
rate function of the distribution at time zero is selected by taking the right 
or left limit of the conditioning. 



and, similarly 



lim lim F(XZ = -LX2 = c) = 5 a 

c->0,fe>0 T->oo 



lim lim P(JCJ = .|A? = c) = «y_ a 

c— s>0,£><0 T— >-oo 



Summarizing our findings, let us denote 3@t the set of bad configurations 
then we have 

Theorem 3.1. 1. Short times: no bad configurations. 
For T < ^ T = 0. 

2. Large times: unique bad configuration. For T > , £$ T = {0} 

3.2 Brownian motion with constant drift 

The case of Brownian motion with constant drift V > is treated similarly. 
The Euler-Lagrange trajectories are once more linear in t, but the cost is 
now 

i(7o) + \ I (is- Vfds 
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which for 74 = A + Bt ending in 77- = b can be computed explicitly and gives 

of which a similar analysis can be given. In particular, choosing b = VT we 
see that VT is the cost is identical to the zero drift case conditioning to be 
at zero at time T, and hence this is a bad point for T > T cr i t , where T cr u 
is the same critical time as for the zero drift case. The analysis around this 
bad point is identical. Notice that the "limiting deterministic dynamics" is 
x — V and the bad point x spec = VT is precisely where this dynamics ends 
up at time T when started from zero. 

3.3 Other rate functions for the initial measure and 
corresponding behavior of Brownian motion 

We now consider other possible scenarios for different rate functions associ- 
ated to the initial measure, and for the Brownian motion with small variance 
as dynamics. The starting measure fi n (dx) satisfies the large deviation prin- 
ciple with rate function i(x). As a consequence, the minimizing trajectory 
to arrive at position b at time T is 74 = Bt + A with B = (b — A)/T and has 
cost 

*b,T(A) = £^¥- + i{A) (15) 
The following scenarios can then occur 

1. i(A) is strictly convex: no bad configurations. Indeed, in that 
case (S'(A) is also strictly convex (as a sum of two strict convex function) 
and hence has a unique minimum. In this scenario, there are no bad 
configurations, and the optimal conditioned trajectory is always unique. 
This corresponds to "high temperature initial measure" and "infinite- 
temperature dynamics" , which always conserves Gibbsianness. 

2. Initial field: loss without recovery, with a "compensating" bad 
configuration. As an example we can take i(A) = (A 2 — a 2 ) 2 + A + r. 
For a > 1, this rate function has one local minimum in the vicinity of 
x = a, a maximum in the vicinity of x = and its (absolute) minimum 
in the vicinity of x = —a. This corresponds to an initial field (favorizing 
the minimizer x = —a). i(A) with a = 2 with r = 2.01539 is plotted in 
Figure [2j The minimization of ^^(A) leads to the equation 

AA 3 + 2a(a,T)A= ^- 1 (16) 
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By an analysis of (16) similar for (14), we obtain that there is no bad 
point when T < T crit = ^j, but b = T is bad for all T > T crit . The 
bad point "compensates" the initial field, and therefore has to become 
larger (and positive) when time T increases. 

3. Non-symmetric rate function. To see that the symmetry of the 
initial rate function is not a necessary requirement to produce bad 
configurations, we have the following example. Let i(A) = 7A 6 — 24A 5 + 
9A 4 + 38A 3 -A2A 2 + A0 (see Figure [3]). This rate function has two global 
minima at x = —1 and x = 2 and one maximum at x = 0. The cost 
function corresponding to trajectories arriving at b at time T is 

g htT (A) = 7 A 6 - 2AA 5 + 9 A 4 + 38A 3 + ( -J- - 42 ) A 2 - %A + %- + 40 



2T ) T 2T 

For fixed b, and T large enough, this function has two local minima, 
located at A l {b,T) < A 2 (b,T). Let us denote, for fixed T, 

D T (b) = ^ T (A\b,T)) - £ b , T {A 2 {b,T)) 

If as a function of b, Dj- changes sign, by continuity, there must be a 
value of b* where Dj-(b*) = 0, i.e., where the minima of <£&* t are a t 
equal height. This b* is then a bad point at time T. For T = 1 we 
have £> T (0.499) ~ -0.00182497 < and £> T (0.4999) ~ 0.000868034 > 
0, so at T = 1, there is a bad point at b* E (0.499,0.4999). We 
observe that b* is T dependent and tends to 0.5 as T increases. From 
numerical oberservations, we have b* G (0.4999, 0.49999) for T = 4, 
b* G (0.49999, 0.499999) for T = 39 and b* G (0.499999, 0.4999999) for 
T = 1000. 

General symmetric rate function. For any rate function i(A) which 
is symmetric with respect to x = and which has minima for A ^ 0, 
b = is bad when T is large enough. Indeed, the cost to arrive at is 
from (15): i(A) + ^ which has a non zero minimum as soon as T is 



large enough. 

5. General short-time Gibbsianness. For every rate function i which 
is twice differentiable and its second derivative is continuous and bounded 
from below, we show that for T small enough there is a unique mini- 
mum Ab of Sbpr{A). This is the analogue of "short-time" Gibbsianness 
obtained in the lattice case via cluster expansions [UJ or conditional 
Dobrushin uniqueness [12J and can be proved as follows. 



11 



We look at the equation (see (15)) 



i'{A) = ~ + ^=:f{A). (17) 



Put d = inf A.i"(A). Then we conclude, for 

-4 



that ( 17) has only one real solution A b . Indeed, look at any two adjacent 



intersection points A\ and A 2 of i'(A) and f(A) if there were more than 



one real solution for (17). By the intermediate value theorem, we get 



mm(i"(A l ),i"(A 2 )) < -i < d = Mi" {A). (19) 

J. -A 

This is a contradiction. And further because i"(Ab) > — ^, we have 

4[ T {A b ) = i"(A b ) + i > 0. (20) 

Therefore A b is a minimum. 

6. Non-Gibbsianness for all times. An example where b = is bad 
for all T > is i(0) = 0, i(A) = f^ 1 \x\ cos 2 \dx for A ^ 0. This follows 
from the facts that i"(A) is not bounded from below when A — > and 
i(A) is symmetric about A = 0. To see that indeed for all T > is 
a bad point, we see that the line f(A) = —A/T always intersects the 
graph of the derivative of the rate function. 



4 Ornstein-Uhlenbeck process 

As a second example, we consider the process X t n to be the solution of 

dX t = -KX t + —^dB t 
y/n 

and the initial point distributed as in the previous section, in ([7]), (]8]). 

The cost function for the large deviation principle of the trajectories now 
becomes 

J?{i) = z( 7o ) + (^ + Klsfds (21) 
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The Euler-Lagrange trajectories extremizing | JqHs + K-ls) 2 ds are given by 

7t = Ae K * + Be~ Kt 

by the terminal condition jt = we have 

7t = -Be^V* + 5e~ Kt 

the cost function for such a trajectory can then explicitly be evaluated and 
gives 

4,r(£) = ci5 4 + c 2 £ 2 + c 3 (22) 



where 



Cl = (l-e" 2K *) 4 



c 2 = (-2a 2 (l - e-™) 2 + Ke~ 2Kt (l - e~ 2Kt )) 

c 3 = a 4 (23) 

A similar analysis as in the previous section can now be started. We have a 
unique minimum at B = of the cost function $ for 

1 , ( 2a 2 \ 

T<T crit :=--log — — 24 
2k \2cr + k J 

and for T > T cri t, becomes the unique bad point for this process. 

The cost of an optimal trajectory ending up at h at time T can also 
be expressed as a function of the starting point 70, which gives the explicit 
expression 



<V(7o) = »(7b) + ^rr(7o - &e KT ) 2 ( 2:, ) 



4.1 Ornstein-Uhlenbeck process with constant exter- 
nal field 

The equation for the process Xj, then reads 

dX! 1 = (-KX t + E)dt+ -f=dB t (26) 

y/n 

where E > is a constant representing a (constant) external field. As rate 
function of the initial measure we choose as before (TsJ) . The cost of the 
trajectory is now given by J* L(-f s , A f s )ds with L( 7s ,7 s ) = (j s + k 7s - E) 2 . 
The Euler-Lagrange trajectories are of the form 

7t = Ae Kt + Be~ Kt + — 

K 
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Figure 4: A limiting process (the purple line) with a = 2, k = 0.7, T = 30, 
hence 7^ ~ 2.0, and a corresponding conditioned process (the blue line) with 
E = 0.1, b « 0.142857 



The trajectory cost of an Euler-Lagrange trajectory is given by 2A 2 (e 2KT — 1). 
From this, we derive that the total cost of a trajectory to end up at time T 
in 77- = b is given, as a function of 70, by 

-(7o-(^-f)e--f) 2 



4t(To) = i(7o) + 



3 2kT 



- 1 



The same analysis can then be performed. The "critical" time at which a 
unique bad point starts to appear is the same as in the zero field case, i.e., 



given by (24). This bad point is given by 

. E , 



-kT\ 



(27) 



which corresponds to the point at which the deterministic evolution x t = 
—KXt + E arrives when starting from xq = 0. Notice that total cost to arrive 
at this bad point b is given by 



i(7o) + 



«7o 



e«" - 1 

which is symmetric around 70 = 0. Moreover, for T large the path cost 

2 

contribution which is equal to e2 ^°_ 1 vanishes exponentially fast, and hence 
for large T two minima exist. 
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The corresponding optimal trajectories to arrive at the bad point b are 
starting from 

7o± = ± \/ a2 ~ 2p^iy 

and explicitly given by 

/, E\ sinhf/ci) / E\ smhUCT - 1)) E 

7* = U - - ■ j 7 ^ + 7o - - . A ^ + - 

\ k J smh(Ki ) \ k J sinn(Ki J k 

The trajectory with plus resp. minus sign can be selected by conditioning to 
arrive at b + > b, resp. b~ < b, and letting b + — > b, resp. b~ — > b. Here we 
plot a limiting process with a = 2, k = 0.7, T = 30, hence 7^ ~ 2.0, and 
a corresponding conditioned process with E = 0.1, hence b ~ 0.142857, see 
Figure [4) 



4.2 General drift. 

Let us now consider the process X™ with a general drift /(a;) and variance 
-, i.e., the solution of 

dX t = -f(X t )dt + ^j=dB t 

We assume / : K — > E to be Lipschitz, and odd: f(—x) = —f(x). For the 
rate function of the initial point Xq we choose as before ([7]), ([8]). The rate 
function of the trajectory is now given by 

^{l) = \j\is + f{ls)fds (28) 

and the minimization problem for the optimal trajectory ending at zero 77- = 
becomes now to find 



argmin{^( 7 ) + i( 7o ) : 7t = 0} (29) 
The Euler-Lagrange equations for miminal cost trajectories are given by 

d 2 7 

These equations correspond to classical motion in a potential U satisfying 
U' = —ff, which gives as a possible choice U = —\f 2 - Notice that this 
formal potential U has no physical meaning, but we need it if we want to 
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translate the framework of the Euler-Lagrange equations to Hamilton equa- 
tions. Indeed, the corresponding Hamiltonian is 

2 

H{p,q) = P --U{q) (30) 
In particular, under the Euler-Lagrange equations, 

| - \{f{lt)f = E (31) 

is a constant of motion. Further, we have the open-end and terminal condi- 
tion 

»'(7o) = 7o + /(7o) 

It = (32) 

We can think of these equations as having 70 and E as parameters. The 
terminal condition gives then a relation between E and 70. Notice that the 
trajectory of zero-energy, E = 0, 7 = is always a solution since /(0) = 0. 
We want to show that under some reasonable assumptions, for T small, it is 
the only solution. For this we make the following assumptions. Call S^t{E) 
the collection of all trajectories 7 : [0, T] — >■ R ending at 0, i.e., with 77- = 
and with "energy" E, i.e., such that 

I - Iw^)) 2 = E 

for all < t < T. We impose now the following conditions. 

1. There exist a function if : R — >■ [0, 00) and T > and a constant C > 
such that ip(0) = 0, <p(E) > for E ^ such that for all T < T and 
for all 7 G ^r(-E), 7o7o < 0, 

| 7o | > <p(E) (33) 

and 

|7o| < C V (E)T (34) 

2. The drift function / is locally monotone around 0, i.e., there exist x 
such that / restricted to [0, x ], [— x ,0] is monotone. 

The first condition states that if T is small, and one wants to end at 77- = 
from 70 > 0, then the derivative at zero should be negative, or vice versa. 
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The second part of the condition states that there exist lower bounds for the 
derivative and upper bounds for 70 . 

Coming back to the previous examples: for the Brownian motion case, 
for all 7 G y T (E) we have j t = ±y/2E(t - T) hence and for 7 G S" T {E) 
we have 70 = ^y/2ET, 70 = ±\/ZE, and we can choose <p(E) = \/2E. For 
the Ornstein-Uhlenbeck case we have 7* = B{eT Kt — e~ K ( 2T ^), E = —2ABa 
and if 7^ = we find 70 = \J2E / K^mh^T), 70 = — ^/2E / k cosh(KT) which 
clearly satisfies the conditions, with the ip = ^2E/k. 

The open-end condition requires 

7o + /(7o) = 47o(7o - ° 2 ) 
Hence, for 7 G ^t^E) such that 70 > 0: 

- (p(E) > 70 

> 4 7o (7 2 - a 2 ) - f(C<p(E)T) 

> AC V (E)T(C 2 V (E) 2 T 2 -a 2 )- f(C<p(E)T) (35) 

which is clearly a contradiction for T sufficiently small. Hence for T suffi- 
ciently small, there do not exist E 7^ with 7 G St(E). As a consequence, 
under these assumptions, for small T the zero trajectory is the only solution 



of the minimization problem (29) 



For large times, if we assume that the drift is such that from any starting 
point one can travel to the origin at arbitrary small cost if one has sufficient 
time, i.e., for all xo > 0, 

lim inf 11 (% + f(>y) s ) 2 ds : 70 = x , 7r = \ = 
Uo J 

then this implies that for T large enough that there exists Xq 7^ and a 
trajectory 7 starting from x such that i(xo) < *(0)/2 and 

(% + f{l)sfds : 70 = x ,7t = oj < z(0)/2 







this trajectory 7 clearly has lower cost than the zero trajectory, and by 
symmetry, —7 is a trajectory with identical cost. Therefore, becomes a 
bad point. 



5 Approximately deterministic walks in d = 1 

An "approximately deterministic random walk" is a continuous-time random 
walk with small increments performed at high rate, i.e., a random walk 
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on R that, starting at Xq = x makes increments of size ±l/iV with rates 
Nb(x), resp. Nd(x). In other words, X^ is a Markov process on R with 
generator 

JZ N f{x) = Nb(x) (f (x+j^- + Nd(x) (f (x - - f(x)j 

(36) 

Such walks arise naturally in the context of population dynamics see e.g. 
[7]. The notation b(x) and d(x) is also reminiscent of this interpretation and 
we will call these quantities birth resp. death rates. 

We ask then the same large deviation question, i.e., we start the process 
X^ from an initial distribution fi^ satisfying the large deviation principle 
with rate function ^ -or some natural modification of it if we have to restrict 
the state space- and look for the minimizing trajectory(ies) that end at time 
T at the origin (or at a more general bad point if the dynamics has a drift 
see later). 

The large deviation function for the trajectories can be computed us- 
ing the Feng-Kurtz scheme, i.e., denoting fp(x) = e Npx we computes the 
Hamiltonian 

JT(x,p) = jim i (±(Sf N ff)^ (x) = (e" - l)b(x) + (e^ - l)d(x) (37) 

and the corresponding Lagrangian 

L(x, v) — sup (pv — Jt?(x,p)) (38) 

For the trajectories of {X^ : < t < T}, we have 

F{X N w 7) w e- N ^ L ^ a) ds (39) 

where the informal notation has to be interpreted as usual in the sense of the 
large deviation principle. 

The equations for the optimal trajectories, i.e. for the minimizers of the 
"action" 

^(7)=«(7o)+ / L( ls ,j s )ds (40) 
Jo 

can now more conveniently be written in terms of the Hamiltonian (the 
Lagrangian is a more complicated expression to deal with). 
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Introducing the canonical coordinates (x,p) we have the Hamilton equa- 
tions, together with the terminal condition and the open-end condition cor- 
responding to the choice of the distribution of Xq . 

dH 

x t = {x t ,Pt) = b(x)e p - d(x)e~ p 
op 

a tt 

Pt = - 1 -(x t , Pt ) = -b'(x)(e p -l)-d'(x)(e- p -l) (41) 



dx 



with conditions 



x T = 

Po = i'(x ) = 4x (xl - a 2 ) (42) 
The total "energy" is a constant of motion along minimizing trajectories, so 



we put H(x,p) = E and we can rewrite the Hamilton equations (41) 



E + b(x) + d(x) + x = 2b(x)u 

E + b(x) + d{x) - x = 2d{x)u- 1 (43) 

where u = e p . Which leads to 

x 2 = E 2 + 2E(b(x) + d{x)) + (b(x) - d{x)) 2 (44) 

So we can think now of the cost of a trajectory as a function of two pa- 
rameters: the starting point and the energy (xq,E). Zero-energy corre- 
spond to the "typical trajectory" following the limiting differential equation 
x = b(x) — d(x), which means that the cost of the Lagrangian part of the rate 
function is zero, and only the cost due to the starting point xq has to be paid. 
Non-zero energy trajectories have a strictly positive cost of the Lagrangian 
part of the rate function. The additional terminal condition Xt = b will 
eliminate one of these variables (e.g. E), so that we can think of the cost of 
the trajectory as a function of the single variable (e.g.) x$. 
We now concentrate on three important particular cases. 



5.1 Constant birth and death rates 

If b and d do not depend on x, then the equation for the momentum shows 
that p t = C, hence we have linear Euler-Lagrange trajectories, and corre- 
spondingly the same analysis and phenomena as in the Brownian motion 
case of the previous section. 
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5.2 Mean- field independent spin flip 

A special case, corresponding to independent spin flip dynamics is b(x) = 
(1 — x), d(x) = (1 + x). Moreover, the x-variable is now restricted to [—1, 1]. 
As in the case x G R we assume that initially, xq is distributed according to 
a measure ^ n {dx) on [—1, 1] satisfying the large deviation principle with the 
non-convex rate function (|8j for x & [— 1, 1] and +oo otherwise. In particular, 

oe(o,i). 

The Hamilton equations then read 

x = -x{e p + e~ p ) + e p - e~ p 

p = e p -e~ p (45) 

Taking the derivative w.r.t. time of the first equation and using the second 
equation leads to elimination of p, and the simple second order equation for 
x: to 

§ = 4, (46) 

with solutions 

x(t) = C ie 2t + C 2 e~ 2t 

where C\,C% are determined by the open-end condition and the terminal 
condition. This case was treated before in the context of the Curie- Weiss 
model subjected to independent spin flips in [TT] . 

The equation for the momentum can be integrated and gives 

tanhfe/2) = ±Ce 2t 

Furthermore, since 

E = (1 - x){e p - 1) + (1 + x){e~ p - 1) 

is a constant of motion, we find as possible solutions for x, using that xt = 0: 

x t = ±y/E/4(l + E/4) (e 2 ^ - e 2 ^) 

In particular, as in the Brownian motion case, the zero-energy trajectory 
(E = 0) yields Xt = 0. The relation between the energy, initial position and 
initial momentum is 



Po = log 



'2 + E+ y/(2 + ff) 2 -4(l 
2(1 -x ) 



Zero-energy thus corresponds to zero initial momentum and zero initial po- 
sition. 
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In general, the initial points are symmetrically distributed around the 
origin and related to the energy via 



±y/E/A{l + E/A) (e" 2T -e 2T ) 



Whether or not a non-zero energy solution can be the minimizer is determined 
by the open-end condition: 

p = i'(x ) = 4x (xo - a 2 ) (47) 

This can be viewed now as an equation for E. For small T > 0, 

x = x (E, T) « C{E)T, po = p {E, T) « cE 



which implies that a non-zero energy solution of (47) can not exist for small 
T. For large T, a non-zero energy solution exists yielding two symmetrically 
solution for xq. 

Alternatively, the trajectory cost Cr(7o) of a trajectory starting at 70 
ending up at time T at b = has the following important properties 

1. Symmetry: C T (-j ) = C T (7o) 

2. Small time behavior: limr_>.o Ct(to) = 00 for all 70 7^ 

3. Large time behavior: limT^oo CV(7o) = for all 70 

From these properties it follows that for small T there are no bad points, and 
for large T zero is the unique bad point. Notice that contrary to the Curie 
Weiss model situation analyzed in [9] there are no non-neutral (non-zero) bad 
configurations due to the fact that the rate function of the initial measure is 
here simply a fourth order polynomial. 

5.3 Independent spin-flips in a field 

This corresponds to the choice b(x) = 7(1 — x), d(x) — (1 + x), x e [—1, 1]. 
Here 7 > 1 corresponds to a bias in the plus direction (positive magnetic 
field). The limiting deterministic trajectory is given by 

^ = -(l + 7)*t + (7-l) 

Xt = xoe-^ + ^— (1 - e-( 1+ ^) (48) 
1+7 

This is the zero-energy trajectory starting from xq. 
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Using (44) we find that for a given energy E, the solution for x is of the 
form 

x t = x(E, C, t) = C ie * (1+7) + C 2 e-* (1+7) + C 3 (49) 



with 



C 2 



C(2 + 2 7 + E) 
( 7 + l)2 
£7 

-C(l + 7 ) 2 ) 



, 3 - (^f^) 

where C is an integration constant. 

Remark 5.1. i. Remark that for E = C 3 = (7 — 1)(1 + 7)^ which 
corresponds to the limiting value of the zero energy trajectory. 

2. If 7 = 1, and E 7^ we find C 3 = and recover the solution of the 
form C\e 2t + C*2e _2i corresponding to the optimal trajectories of the 
independent spin flip dynamics. 

The general form of an optimal trajectory arriving at time T at xt = b 
and starting from Xo = 70 is 

„ . sinh(5t) ^ sinh(5(T — t)) _ 

smh(oi ) smh(oi ) 



with 5 = (1 + 7) and where C3 is given in (50). Notice the analogy with the 
case of the Ornstein-Uhlenbeck process in a constant field (27). As in that 
case, the bad point is time-dependent and given by 

b = 1^1 (1 - e~ ST ) 
7 + 1 v ; 

which is the point at which the limiting deterministic dynamics arrives at 
time T when started from xq = 0. The trajectory cost CV(7o) to arrive at this 
bad point satisfies the same properties as the trajectory cost CV(7o) of the 
previous subsection (zero field case). Hence, for T large two minimizing 70 of 
the total cost function appear which correspond to two optimal trajectories. 
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